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Abstract 

This paper solves the dynamic portfolio choice problem. Using an explicit solution with a power 
utility, we construct a bridge between a continuous and discrete VAR model to assess portfolio sen¬ 
sitivities. We find, from a well analyzed example that the optimal allocation to stocks is particularly 
sensitive to Sharpe ratio. Our quantitative analysis highlights that this sensitivity increases when 
the risk aversion decreases and/or when the time horizon increases. This finding explains the low 
accuracy of discrete numerical methods especially along the tails of the unconditional distribution 
of the state variable. 
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Numerical solution. 

JEL Classification: Gll; G12. 


Introduction 

Since at least Merton (1971), many results on portfolio optimization problems have been obtained 
in a continuous time framework. It is still difficult to solve optimal portfolio problems when there 
is some degree of predictability in asset returns, i.e. when the investment opportunities are time- 
varying. A great number of papers have proposed to use a VAR model to forecast returns and study 
its implication for the long-term portfolio choice problem. As a result the academic literature has 
followed two main directions. The first one relies on mathematical tools and establishes some explicit 
solutions (see among others Kim and Omberg (1996), Liu (2007) and references therein). The second 
line of research consists to implement some challenging numerical methods. In fact, Barberis (2000) 
developed a discretization state space method which serves as a benchmark. Brandt et al. (2005), van 
Binsbergen and Brandt (2007), Garlappi and Skoulakis (2009) among others use some sophisticated 
backward induction techniques and evaluate the accuracy of their results by comparing them to the 
discretization state space benchmark. However, all discrete numerical procedures approximate directly 
or indirectly a highly non linear value function and cannot explicitly separate the so-called hedging 
demand from the so-called myopic demand. Garlappi and Skoulakis (2011) provide a general discussion 
about approximations accuracy in discrete time. 
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AErudite, UniversitAl’ Paris-Est, D j ibril. Togola@u-pec. f r. 
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This paper works in continuous time and uses the explicit solution for the portfolio choice problem, 
then constructs a bridge between continuous and discrete VAR model as in Campbell et al. (2004). 
In fact, these authors provided evidence that there should exist minor discrepancies between results 
under discrete and continuous time models. Thus, numerical results that we derive from continuous 
time are indirectly comparable to those of Garlappi and Skoulakis (2009). We show that, for large 
degrees of risk aversion and/or small horizon, when the state variable closes to its unconditional mean, 
the two numerical results are quite similar. Otherwise, results under our explicit solution in continuous 
time exhibit some discrepancies with Garlappi and Skoulakis (2009) when the risk aversion decreases 
and/or when the time horizon increases. We argue that this is due to large sensitivity of total demand 
to state variable (Sharpe ratio). 

The paper is organized as follows. Section 1 exposes the way we map the continuous-time invest¬ 
ment opportunity set and the discrete-time one. Section 2 gives some insights on the explicit solu¬ 
tion for the long-term investor with CRRA preferences. Section 3 gives some numerical results based 
on Brandt et al. (2005) example. 


1 Investment opportunity sets 

We first expose the model in a continuous-time framework and in a discrete-time framework to study 
the impact of a predictable component in stock returns. Next, we show how to recover continuous-time 
parameters that are consistent with discrete-time VAR estimates. 


1.1 Opportunity set in continuous time 


We start by assuming that two assets are available for the investor (Gampbell et al. (2004) and Kim and 
Omberg (1996) among others). On the one hand, the riskless asset pays back a constant interest rate r 

pf 

-j = rt (1) 

Pt 


where p( denotes the price of this asset at time t. On the other hand, there is a risky asset whose 
price Pf satisfies the following diffusion process 


ft 

Pt 


^ ft 




( 2 ) 


where denotes a scalar Brownian motion with zero drift and unit variance rate. The drift rate ^ 
follows a diffusion process as well. It is supposed to be time-varying and state variable dependent. The 
volatility of the risky asset is assumed to be constant. This is not a strong assumption for the long-term 
investor (see Gampbell and Viceira (2002)). Let denote the Sharpe ratio i.e. the market price of 
risk/reward for buying/selling one unit of risky asset 



(3) 


Then the Sharpe ratio is assumed to follow the usual “Ornstein-Uhenbeck” diffusion process 

X, = (-X,)t + B^ ,,>0 (4) 

where B/ denotes another scalar Brownian motion with zero drift and unit variance rate. Parameters 
and denote respectively the unconditional mean and the mean reverting parameter of the Sharpe 
ratio Xf.. In fact, parameter reflects the rate by which the shocks on Sharpe ratio dissipate and then 
reverts towards its long-term mean . Finally, parameter denotes the instantaneous volatility of Sharpe 
ratio. It controls the diffusion rate of the process. 
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The above equations imply that instantaneous return on stocks Pj/Pf follows a diffusion process 
whose drift is mean-reverting and whose innovations are correlated with those of the market price of 
risk itself, with the correlation coefficient. Thus the following equations hold. 

Pf/Pf = (X, + r)t + Bf 

with = t. Equations (5) and (6) define a joint stochastic process in continuous time 

1.2 Opportunity set in discrete time 

The standard model in discrete time is a restricted VAR(l) process which captures predictability of 
stocks returns (see Barberis (2000) for instance). We focus on the example analyzed in Brandt et al. 
(2005) that was reused in van Binsbergen and Brandt (2007) and in Garlappi and Skoulakis (2009). 
The log excess returns of the risky asset InP^+i — are assumed to be predictable by the log dividend- 
to-price ratio (r^ denotes the risk-free rate and is equal to 6% in annualized basis). The joint 
dynamics of these two variables are specified such that 

InPf+i-r-^ ^a, + b,Zc + [^^ 

Zt+i + + t+i 

with 

1+1 I ~iV 

t+ij 

In fact, Campbell and Shiller (1988) forcefully claim that, if returns are predictable, at least, the log 
dividend-to-price ratio should capture some of that predictability. A substantial long-standing empirical 
literature has documented many properties of these two regressions. Brandt et al. (2005) report the 
following estimated values (using the CRSP U.S. quarterly index from January 1986 to December 1995) 

b, = 0.060 b^ = 0.958 — - -0.941 

r z 

The returns are weakly predictable, the dividend yield is highly persistent and the shocks are strongly 
negatively related. 



(7) 

( 8 ) 


(5) 

( 6 ) 


1.3 Recovering continuous-time parameters from discrete-time VAR 


We closely follow Campbell et al. (2004) to recover the parameters of the continuous-time system 
eqs (5)-(6) from the restricted VAR(l) eqs (7)-(8). However, Campbell et al. (2004) use the risk 
premium as state variable; we prefer to use the Sharpe ratio. In matrix form, the discrete-time VAR 
eqs (7)-(8) is 


InPf+i 

2 t+l j 




( 10 ) 


The first step is to aggregate the continuous-time model over a span of time taking point obser¬ 
vations at evenly spaced points with t = We then obtain, using the 

discretization method developed by Bergstrom (1984) 



(_2/2 + )t_(i_-t)A 

(1--^) J 



( 11 ) 
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Discrete-time world 

Continuous-time world 

Brandt et al. (2005) 
InPt+i-rl = a, + b,Zt + [^^ 
Zt+i =az + hzt + ?+i 
? 

v(?) = .^ 

Cov(^4) = r. 

Models 

Kim and Omberg (1996) 

p{ /P{ = rt 

PjP, = iXt + r)t + BP 

A, = (-Ajt + Bf 

BfB^ = t 


Parameter values 

Brandt et al. (2005) 

Our computations eqs (13)-(18) 

rf 

0.015 


a. 

0.227 

r 0.015 

K 

0.060 

0.111 

az 

-0.155 

0.0429 

bz 

0.958 

0.0060 

2 

r 

0.0060 

0.0542 

2 

z 

0.0049 

-0.941 

rz 

-0.0051 



Table 1: Recovering continuous-time parameters 


where 



( 12 ) 


with -1- \l 1—^Z^ where Bf and are two independent brownian terms. 

The second step is to apply a linear transformation for the process in (11) so that we can relate 
the parameters of the transformed system to the parameters of the matrix form (10) of the discrete¬ 
time VAR model. Thus, when we normalize time span t = 1, since everything is in quarter, we get 
(for b^, > 0) 


rf 

(13) 

a^b^ a^ + ^^/2 

r(l-M^ r 

(14) 

-ln(bj 

(15) 

r 

(16) 

K- 

(17) 

r 

rz 

(18) 


r Z 


The appendix proves these results. Table 1 shows the value of the parameters of the continuous-time 
equivalent VAR implied by the Brandt et al. (2005) estimates. 

2 Portfolio choice problem in continuous time with CRRA preferences 

We can now solve the portfolio choice problem of the investor with long-term horizons who faces to the 
investment opportunity set described in the previous section. We rely on the recent advances in Honda 
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and Kamimura (2011) who use the verification theorem and show that the explicit solution provided 
under continuous time is in fact an optimal solution especially for risk aversion greater that one. 

We consider an investor with initial wealth > 0 who has only two assets (riskless short-term 
bond and stocks) available for investment. The financial markets are incomplete. Furthermore, the 
investor can undertake continuous trading, he has no labor income and only cares about terminal 
wealth Wj where T is the finite planning horizon. The dynamics of price changes are described by (1) 
and (5)-(6). If ^ is the share of wealth invested in stocks, the instantaneous wealth would be given by 


Wt 




(19) 


Properly substituting the dynamics of Pf/Pj and /P^, wealth dynamics (also called the budget con¬ 
straint) becomes: 

Wt-(tXt-Fr)Wft-FfWtPf (20) 

Notice that wealth process reflects uncertainty in instantaneous returns (term Pf) and about the state 
variable (the term Xf.). Given this formalization about wealth process, at time to, the investor’s opti¬ 
mization problem can then be expressed as 


max ^u(W-j-) subject to the constraint (20) fixed (21) 

to “ “ 


where Ej denotes the operator of conditional rational expectation at date tg, the time discount param¬ 
eter (with > 0) and u() the utility function defined over terminal wealth. Let tg) defines 

the value of the problem defined in (22) at time tg 

tg) - max E,^(22) 

to 


The Bellman equation generalizes this problem for every time t so that 


J(Wt,Xf, t) - max E,J{W^+Wt,X^+X^, t+t) 


(23) 


Equation (23) emphasizes the fact that current optimal decisions depend on the conditional expected 
value of the problem which, in turn, is intimately linked to future wealth and the state variable. Apply¬ 
ing Ito’s lemma to the Bellman equation, we find that 


0 = max 


4-iitXt + r)W, + ‘^ + ^(-X,) + 


lW, 


1 V 


1 V 


-I _ 

2 ^ ' 2 


^ + - tWt 

9 T7- JAr -w ^ f 


t ^ 


(24) 


The first order condition of equation (24) with respect to ^ implies that 


^ _ J/W, 1 X, V/(W,X,) 1 


(25) 


Merton (1971) was the first to propose such additive decomposition between a myopic demand (first 
term) and a hedging demand (second term) of the optimal allocation to stocks. There is no hedging 
demand especially when the opportunity set is nonstochastic ( = 0) or when the opportunity set is 
uncorrelated with asset returns ( = 0). 

Now, we need to explicitly define the function J(). The first conjecture (see Kim and Omberg 
(1996)) is to assume 

J(1V„A„ t) - -'u(lV,) [/(A„ t)] (26) 
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where /() is an auxiliary function with the terminal condition f{Xj-,T) = 1. We consider the CRRA 
preferences uiW^) = wl~/{I—') where is the coefficient of relative risk aversion. Thus, the hedging 
demand in (25) could straightforward be expressed as 

fl^t _ In/ 

f ~ X, 


Then, under CRRA preferences hypothesis, the optimal allocation to stocks could be expressed as 

In/ 


So, the Bellman equation (24) can be rewritten as 

.2 


(27) 


iL 

7 

/, 


1 - n fii- 


0=7 + —'-- +^x^ + -T—Xt + -f(-Xt) + 


f 


f 

_ ' 

/ 


// 2 


' f 

j V 


—-+ — —n- 1 ) 
/ 2 I / i 2 ^ 


(28) 


where we use intuitive notations for the derivatives of the function /(). Equation (28) is a partial differ¬ 
ential equation which admits analytical solutions especially if utility is logarithmic ( = 1 by FHopital’s 
rule) or if markets are complete ( = ± 1 ). 

The second conjecture is to assume 


/(Af,t) = exp 


Co(t) + Ci(t)X, + ic2(t)X 


2 

t 


(29) 


where Co(t), Ci(t) and C2(t) are some undetermined time varying coefficients (with Co(T) = Ci{T) = 
= 0). Under this conjecture, using equation (27), the optimal allocation to stocks is 


★ 

t 


IX. 

-^ + [Ci(t) + C 2 (t)X,] 


(30) 


We only need to recover Ci(t) and C 2 (t) coefficients. 

This conjecture was also used by Kim and Omberg (1996) and by Liu (2007) among others. More 
recently, Honda and Kamimura (2011) show that the explicit solution derived from the Bellman equa¬ 
tion is in fact, even if the markets are incomplete, an optimal solution to the problem of the long-term 
investor who only care about terminal wealth and who have a risk aversion larger than unity. 

Let us substitute our second conjecture in the equation (28) 



62 2 


[Cl b 1 

0 - 

- h d C2 b Cn + C 

t 2 

Xl + 

- 1 — Cl -l- C 2 + fl C1C2 

t 2 


Co 1 - ^ a 2 

-1-r — -|- C;^ H— C 2 H—C, 

t ^ 22 ^ 


(31) 


where a = [1 -I- (1—)(^—1)] b = 2[(1—)/ — ] and c = (1—)/^. As, whatever the value of X^, the 
equation (31) must hold, all terms within brackets are simultaneously set to zero to solve the equation 
for CqO, CiO, and C 2 (). Defining the discriminant D 

D = — 4ac 


one can check that if > 1 then D > 0. Thus, the two solutions of interest are given by 

2 c(l-“(^“^^) 

2 -(b + )(l 
4c (l--(^-0/^)^ 

2 -(b + )(l 


C 2 (t) 

Ci(t) 


(32) 

(33) 
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Figure 1: Myopic (dashed line) and hedging (solid line) demands as function of risk aversion for = 

where denotes -/D. Kim and Omberg (1996) called this the normal solution and discussed about some 
alternative solutions those are beyond the scope of this paper. The appendix provides details about (32) 
and (33). It is easy to see that there is a linear relation between QO and C 2 (). Then, for > 1, Cj 
and C 2 are always negative. As a result, since < 0, the hedging demand is always positive when the 
preferences are not logarithmic (more precisely for >1) and the market price of risk is positive. 

3 Numerical results 

As mentioned above, we illustrate our approach using the well documented Brandt et al. (2005) exam¬ 
ple. Table 1 collects the continuous-time parameters recovered from this example. For comparison 
purposes, we also use the Garlappi and Skoulakis (2009) results, obtained from the same discrete-time 
VAR(l) estimates and by means of a sophisticated numerical method. 

Figure 1 and table 2 help to understand the long-term investor problem. For =1 i.e. the case of 
logarithmic utility, no hedging demand is required. For this case, the dynamic portfolio choice reduces 
to static one whatever the time horizon. Otherwise, for > 1 and horizon longer than one, under CRRA 
preferences and mean reverting returns, agent should have a positive hedging demand to prevent 
adverse changes in investment opportunities (Merton, 1971). However, for —> 00 , more specifically 
for a very conservative agent, stocks are not attractive. Thus, he would not invest into stocks since the 
total demand (sum of myopic demand and hedging demand) converges toward zero. Our results reset 
all theses basic important features. 

The total demand is sensitive to risk aversion. Results from previous studies imply that myopic 
and hedging demands are more sensitive to small values of risk aversion. We confirm this and argue 
that the sensitivity of hedging demand to state variable is maximal near the critical point 2. Our 
equation (30) and figure 1 show this evidence. To quantitatively see this, just evaluate the derivative 
of * with respect to X. 

Table 2 reports the evidence that both myopic and hedging demand are sensitive to initial value of 
Sharpe ratio. These two components of optimal allocation individually increase with the percentile of 
the Sharpe ratio unconditional distribution. Thus, the total demand exhibits the same behavior. This 
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T 


= 5 





= 15 





^(10) 

^(30) 

^(50) 

o' 

^(90) 

^(10) 

^(30) 

^(50) 

^(70) 

^(90) 

10 

MD 

-34.0 

3.0 

28.6 

54.2 

91.1 

-11.3 

1.0 

9.5 

18.1 

30.4 


HD 

-10.9 

3.5 

13.5 

23.6 

38.0 

-4.6 

1.5 

5.7 

9.9 

15.9 

20 

MD 

-34.0 

3.0 

28.6 

54.2 

91.1 

-11.3 

1.0 

9.5 

18.1 

30.4 


HD 

-15.9 

10.8 

29.2 

47.7 

74.3 

-7.2 

4.9 

13.3 

21.6 

33.7 

30 

MD 

-34.0 

3.0 

28.6 

54.2 

91.1 

-11.3 

1.0 

9.5 

18.1 

30.4 


HD 

-16.0 

19.8 

44.7 

69.5 

105.3 

-7.7 

9.8 

21.9 

34.1 

51.6 

40 

MD 

-34.0 

3.0 

28.6 

54.2 

91.1 

-11.3 

1.0 

9.5 

18.1 

30.4 


HD 

-13.2 

29.1 

58.3 

87.6 

129.8 

-6.5 

15.5 

30.7 

46.0 

68.0 


For each risk aversion , the first line reports the myopic demand (MD) and the second line the hedging 
demand (HD) tvithout short selling constraints. We present the results for 5 different initial values 
of the Sharpe ratio X. Each value corresponds to the p-th percentile of the unconditional distribution 
of X, defined by equations (49) and denoted by X^p^ where p takes values 10, 30, 50, 70, and 90 
(thenX(5o) = ). 

Table 2: Myopic and hedging demands for investment horizon of T quarters (%) 


is consistent with Campbell et al. (2004) among others. In fact, high Sharpe ratio or equivalently high 
risk premium relative to volatility signals better investment opportunities. Therefore, optimal fraction 
to allocate into stocks should increase from the knowledge of mean reverting parameter that serves to 
quantify the expected Sharpe ratio. 

Myopic demand is independent from time horizon while hedging demand increases nonlinearly 
with time horizon. However, table 2 quantitative figures suggest that this relation is almost linear but 
small changes in horizon induce substantial hedging demand. Horizon effect is important but quiet 
monotonic for a given percentile of the state variable unconditional distribution. All changes in total 
demand for fixed risk aversion and state variable are due to changes in horizon and are large for small 
risk aversion. 

The horizon effect on hedging demand is important in optimal allocation because it widely domi¬ 
nates for longer horizons. In fact, when horizon is greater than 20 quarters, hedging demand becomes 
always greater than myopic demand when the Sharpe ratio initial value is between 30 and 70 per¬ 
centiles. 

We finally use the common assumption of no-borrowing and short-sale constraints. Thus, in table 3, 
we restrict all portfolio weights between 0 and 1. One can notice that we generally obtain values 
fairly close to those of Garlappi and Skoulakis (2009) while frameworks are not the same. Garlappi 
and Skoulakis (2009) worked in discrete-time and initial values of their state variable are drawn for 
the unconditional distribution of quarterly dividend price ratio. They use a sophisticated numerical 
optimization technique. We work in continuous time (no numerical optimization) and our initial values 
are computed using the unconditional distribution of continuous Sharpe ratio that we discretized in 
points observation and recovered using the same quarterly dividend price ratio. However, a closer 
inspection of table 3 figures show that the optimal allocation to stocks is more sensitive to the state 
variable and to the time horizon than the sensitivity obtained by Garlappi and Skoulakis (2009). We 
run some numerical simulations, within the discrete-time framework, to evaluate our results in order 
to find the causes of the discrepancies between the two frameworks. We were unable to qualitatively 
and quantitatively invalidate our results. 

To test our results, we run some forward pure simulations in discrete time. More precisely, for 
instance, we explore the case where the initial value of the Sharpe ratio is the 30-th percentile (Aq = 
A( 3 o), the relative risk aversion is equal to 5 ( = 5), and the planning horizon is equal to 10 (T = 10 
quarters). With this configuration, when we get an initial optimal allocation to stocks of 0.065, Garlappi 
and Skoulakis (2009) obtain twice as many (0.133, see table 3). That’s large. Thus, we first build a 
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Horizon in quarters (T parameter) 


Figure 2: Optimal allocation to stocks as function of the horizon for = 5 (first panel) and for =15 
(second panel) for 5 different initial values of the Sharpe ratio X (as in table 2 or 3) 



Figure 3: Path of optimal allocation to stocks for = 5, Xq = Xpg), and T — 10 obtained by explicit 
solution (solid line) and by simulation and trial {0.05,0.10,0.15,0.20,0.25} grid (dashed line) 
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= 5 





= 15 





T 


^(10) 

^(30) 

^(50) 

^(70) 

^(90) 

^(10) 

^(30) 

^(50) 

^(70) 

-^(90) 

10 

LT 

0.0 

6.5 

42.1 

77.7 

100.0 

0.0 

2.5 

15.2 

27.9 

46.3 


GS 

0.0 

13.3 

43.2 

73.1 

100.0 

0.0 

4.3 

15.4 

27.0 

44.7 



0.0 

- 6.8 

- 1.1 

4.6 

0.0 

0.0 

- 1.8 

- 0.2 

0.9 

1.6 

20 

LT 

0.0 

13.7 

57.8 

100.0 

100.0 

0.0 

5.9 

22.8 

39.7 

64.1 


GS 

0.0 

24.4 

57.2 

89.7 

100.0 

0.0 

10.7 

25.1 

40.4 

63.2 



0.0 

- 10.7 

0.6 

10.3 

0.0 

0.0 

- 4.8 

- 2.3 

- 0.7 

0.9 

30 

LT 

0.0 

22.8 

73.2 

100.0 

100.0 

0.0 

10.8 

31.5 

52.1 

81.9 


GS 

0.0 

32.8 

68.4 

100.0 

100.0 

0.0 

17.5 

35.2 

54.0 

80.7 



0.0 

- 10.0 

4.8 

0.0 

0.0 

0.0 

- 6.7 

- 3.7 

- 1.9 

1.2 

40 

LT 

0.0 

32.0 

86.9 

100.0 

100.0 

0.0 

16.5 

40.2 

64.0 

98.3 


GS 

0.0 

38.8 

77.6 

100.0 

100.0 

0.0 

24.1 

44.5 

65.7 

94.6 



0.0 

- 6.8 

9.3 

0.0 

0.0 

0.0 

- 7.6 

- 4.3 

- 1.7 

3.7 


For each risk aversion , the first line reports our results (LT aA§ optimal allocation to stocks in continuous 
time), the second line the Garlappi and Skoulakis (2009) results (GS aA§ optimal allocation to stocks in 
discrete time), and the third line reports the difference between our results and Garlappi and Skoulakis 
(2009) results. We present the results for 5 different initial values of the Sharpe ratio X calibrated using 
the same estimates involving dividend price ratio as in GS. Each value corresponds to the p-th percentile 
of the unconditional distribution of X, defined by equation (49) and denoted by X^p^, where p takes 
values 10, 30, 50, 70, and 90. 

Table 3: Optimal allocation to stocks for investment horizon of T quarters (%) 


sample of size 100000 for Zf_i_2, andzj+io and for luPf+i, lnPf_,_2> ^nd InP j+iq using the 
restricted VAR(l) eqs (7)-(8). We choose the grid G — {0.05,0.10,0.15,0.20,0.25} for trial allocations 
to stocks, to overlay both our and Garlappi and Skoulakis (2009) solutions. Then, for each path in 
the sample, the value of terminal wealth is computed from the cartesian product G x G x x G of all 
possible strategies. The computational burden is very high as we evaluate 5^° = 9 765 625 strategies. 
Figure 3 shows that the forward path in discrete time (no numerical optimization) closes to the path of 
our explicit solution particularly at the critical starting point, the 30-th percentile of the state variable 
for small risk aversion ( = 5). 


Conclusion 

We examine the “continuous-time detour” to solve the long-term investor problem when the stock 
returns are predictable. We obtain an explicit optimal solution in the continuous-time world and, 
after recovering the continuous-time parameters from the discrete-time world estimates, we reuse such 
solution to assess the sensitivities of optimal allocation to the initial values of the state variable, to 
the risk aversion and to the time horizon. We find greater sensitivities than those reported in the 
literature. We also find that the sensitivity of total demand to the state variable is not uniform along 
the unconditional distribution of the state variable. 

Previous numerical approximation techniques that deal with the problem we consider are subject 
to some numerical errors. Therefore, they do not always provide accurate results. We show that 
the hedging demand part of allocation dominates at longer horizons and it is very sensitive to state 
variable especially when risk aversion decreases and/or the time horizon increases. This finding could 
explain the low accuracy of discrete numerical methods especially along the tails of the unconditional 
distribution of the state variable. 
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Appendix 

A. 1 Proof of continuous VAR recovering by discrete VAR 

The matrix (10) could be rewritten as 

InPt+t -r^ = a,.-\- b,.z^ -+- (34) 

Zf+t = + f+f (35) 
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Equations (34) and (35) describe a joint process of an econometric model in which z denotes the 
dividend price ratio. The corresponding discretized version of the continuous time model in matrix (11) 
could be rewritten as 


logP - rt = (-2/2 + )t - (1 - -^) + (1 - + UP 


Comparing the expectations of (34) and (36), we get 

Iterating forward (38), t periods ahead and using (35), we obtain 

+ (- 2/2 + )f - (1 - -^ 7 - + (1 -= 

bj. bj. b^ bj. 

““ + (-(-"/2 + )^ + (1 - - (1 - + U 

After some algebra, we find that 


X = 


a^b. 


+ (l~^z)^—^ + (-/2 + )t- 


1- 


+ h„Xf — 


1 - 


-t t+t 


(36) 

(37) 


(38) 


(39) 


(40) 


Notice that limj^o (1 “ ^) — 1- Finally, comparing equation (40) to (37), equations (13)-(15) directly 
follow. To compute the associated second moments, one can compute the variance of U vector in (12). 


Var 



rt 

FF' 

J=o 



(41) 


where 


F = 


1 (1--^) 
0 


0 " 


(42) 


Using a matching procedure involving equations (34)-(41) we can directly reset equations (16)-(18). 
Furthermore, the resulting general formula for every t become such that 


Var(X) =-(l - -2f) 


Cov(X, logPt+f)=-(l- f) + 2 (l- ^0 


(43) 

(44) 


22 ' 


Var (ln(Pt+t)) = |^2 + 2p-+ — J t - 2p- (l - 

22 


22 




(45) 


Where the instantaneous standard deviation of X denoted is given by equation (17). Again, notice 
that, for small , i.e. when t —> 0, the term (1 — “^) —> t. So when t = 1, all second moments 
could be approximated by their instantaneous counterparts. Otherwise, when t ^ 1, these kinds of 
approximations become no longer valid. Campbell et al. (2004, p. 2208) discuss about pitfalls resulting 
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for this case. Taking this into account, for instance, one can compute the terminal conditional variances 
by just setting t = T and t = 0. 

The unconditional moments ofX that have been used in this paper are derived from equation (38) 
when t is normalized to one. 

(46) 

(47) 

Then, the unconditional mean of X is 


a„ + 

Xt =x + -—— 

2 

a,+ b,E(zt)_ a, + b,ajil-bj 

J —— H-— — H- 

^” 2 2 


a^b, ar + r/2 

r 


(48) 


In fact, we have used the result = ^ in equation (16) and the fact that z follows an AR(1) process 
(Brandt et al. (2005) followed by Garlappi and Skoulakis (2009) among others). Thus its uncondi¬ 
tional moments are knovra, E (zj) = a^l{l — b^) and Var (Zf) = ^/(l — bf) (Hamilton (1994, p. 53)). 
So, under equation (38), one can match all unconditional percentiles Zj-^) with their unconditional 
counterparts (p denotes the p-th percentile) and get the optimal policies for those values. Since 
the hedging demand is very sensitive to state variable, we directly draw from the unconditional 
distribution of the point observation Xf. of our continuous time state variable X in order to avoid 
computational errors. As a result, Aij-so) = and the following unconditional distributions hold. 


z ~ Ai 


CL^ 

l-b.’ 


A ~ A 





(49) 


A. 2 Proof of parameters and C 2 

Regarding (31), the solution for C2 could be derived from the following equation: 

C2 , 

-h d C9 b Co + c — 0 

t 2 

This equation could straightforward be rewritten as 

1 

C2 - -(T - t) 


(50) 


CL Cq + b C 2 c 


(51) 


Since parameters a, b and c are constant, given C2(T) = 0, the integral table provides the solution 
for € 2 - Substitute this into the following equation 


Cl 

-h Co -l- I — -l- d Co 1 Cl = 0 


(52) 


that we derived from equation (31). Again using the terminal condition Ci(T) = 0 and the constant 
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variation method, one can get the solution for as follows : 


Ci(t) 


rT 


t 

rT 


t 

rT 




(-C,(s))s 


, ( + 6)((^-)-l) + 2 

2c(r-t)/2 rT 


( + bX(^-f)-l) + 2 
4c/ 

-( + bXl-(r-f)) + 2 

2 2c 

~2-( + bX^^“'^-l) 




( 


(r-0/2 


4^ 


(^1 _(r-t)/2y 


2(1 


. (r-t)/2 


) 


1 _(r-t) 


C2(t) 
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